Strongly inhibited transport of a ID Bose gas in a lattice 
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We report the observation of strongly damped dipole oscillations of a quantum degenerate ID 
atomic Bose gas in a combined harmonic and optical lattice potential. Damping is significant for 
very shallow axial lattices (0.25 photon recoil energies), and increases dramatically with increasing 
lattice depth, such that the gas becomes nearly immobile for times an order of magnitude longer than 
the single-particle tunneling time. Surprisingly, we see no broadening of the atomic quasimomentum 
distribution after damped motion. Recent theoretical work suggests that quantum fluctuations can 
strongly damp dipole oscillations of ID atomic Bose gas, providing a possible explanation for our 
observations. 
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The ability of highly degenerate quantum systems to 
sustain dissipationless flow is one of the most striking 
manifestations of quantum mechanics. However, trans- 
port in such systems can be dramatically modified by the 
presence of a relatively weak, but rapidly spatially vary- 
ing ("corrugated") potential along the transport axis. 
For example, the periodic potential of an optical lat- 
tice inhibits transport in a degenerate Fermi atomic gas 
[H El hut not, in general, in a degenerate Bose gas 
(i.e., Bose-Einstein Condensate (BEC)) However, 
under certain conditions, htehly dissipative transport in 
a BEC in an optical lattice la LZtlj l can arise from nonlin- 
ear dynamical instabilities |9UlOUll|. In low dimensional 
systems, of which ID atomic gases 0, 0, 0, 0] 
and superconducting nanowires are important exper- 
imentally realized examples, a corrugated potential can 
cause dramatic changes in ground state and transport 
properties. 

We study inhibited transport in a ID Bose gas in the 
presence of an optical lattice along the ID axis. In the ab- 
sence of such a lattice, dipole oscillations are undamped 
[T3 |. since it is a general result that the dipole mode of 
a harmonically confined gas is unaffected by 2-body in- 
teractions (generalized Kohn's theorem) ^^|. This result 
does not strictly hold for a combined harmonic and peri- 
odic potential; nevertheless, undamped oscillations have 
been observed in 3D BECs for small amplitudes and weak 
interactions^ . 

In this Letter we report a study of strongly damped 
dipole oscillations of a ID Bose gas in a combined har- 
monic and periodic potential, under conditions for which 
undamped motion has been observed previously for 3D 
BECs. This striking difference between ID and 3D was 
recently reported, qualitatively, in Ref. .13] . Here we 
measure the damped motion as a function of axial lattice 
depth. Significant damping is induced by very shallow 
lattices, and in deeper lattices the motion is overdamped 
to the degree that the gas is nearly immobile for times an 
order of magnitude longer than the single-particle tunnel- 



ing time. We emphasize, and discuss further below, that 
the inhibited transport is not due to Bloch oscillations 
0, |2(| , where transport is frustrated by Bragg reflection 
at the Brillouin zone (BZ) boundary, as has been seen in 
previous experiments 0, Q 0] • 

Our method to realize an ensemble of independent ID 
Bose gases is similar to earlier work ,12). We produce a 
nearly pure 87 Rb condensate of N = (0.8 — 1.6) x 10 5 
atoms in the \F = 1,itif = —1) state in a Ioffe-Pritchard 
magnetic trap [y x = v z = 29 Hz, v y = 8 Hz). We next 
partition the BEC into an array of independent, verti- 
cal ID "tubes" by adiabatically applying a transverse (in 
the xy plane) 2D confining lattice (Ll Ha Hi HI Hrj l2l| . 
The confining lattice is ramped on during 200 ms to a 
depth of approximately 30 Er (where Er = h 2 /2m\ 2 is 
the photon recoil energy, and A is the laser wavelength). 
The combined magnetic and optical potential results in 
approximately 5000 occupied tubes, each with an axial 
frequency of ujq/2tt rs 60 Hz. We observe a Thomas- 
Fermi density envelope in the combined magnetic and 
optical potential, and calculate 0] cloud radii of r x — 
14(1) fxm H3, r v = 20(1) fxm, and r z = 10.6(5) (im for 
N = 1.4 x 10 5 . From this we estimate a peak ID density 
of 4.8(4) x 10 4 cm -1 in the central tube, and a peak 3D 
density of 4.7(4) x 10 14 cm~ 3 . Subsequently, we corru- 
gate the tubes by adiabatically applying, over 20 ms, an 
axial (vertically along z) ID lattice. The Rayleigh length 
of the axial lattice beams is large enough that they do 
not significantly modify the axial harmonic potential. All 
lattice beams derive from a single TbSapphire laser op- 
erating at A = 810 nm, far detuned from the atomic 
resonances at 780 nm and 795 nm. The pairs of lattice 
beams are detuned from each other by 6 MHz, making 
them effectively independent j2^|. The final configura- 
tion consists of three independent standing waves, each 
formed from a pair of counter-propagating beams. 

We excite dipole oscillations of the center-of-mass of 
the atoms in all the tubes by suddenly (< 150/is) ap- 
plying a linear magnetic field gradient, thus displacing 
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FIG. 1: Damped oscillations of a ID Bose gas in an optical 
lattice. Shown are plots of velocity versus wait time t w from 
t w = to 110 ms, and for axial lattice depths of (a) OEr, (b) 
0.25E R , (c) 0.50E r , and (d) 2.0E R , where E R is the photon 
recoil energy (see text). 



the total harmonic trap (but not the lattice) axially by 
zq w Sfim. This displacement is less than 30% of r z , 
and corresponds to approximately eight axial lattice sites 
[23 | . The waists of the Gaussian transverse lattice beams 
(wq « 210 /im) are much larger than both zq and the size 
of the trapped cloud. 

The oscillation in the position of the atoms is too small 
for our imaging system to clearly resolve. We therefore 
observe oscillation in velocity by waiting a variable time 
t w after the initial displacement, then suddenly turn- 
ing off all trapping potentials (with time constants of 
~ 250 /is and w 150/iS for the optical and magnetic 
potentials, respectively), and imaging the atoms after a 
time-of-flight ixOF = 18.4 ms. The turn-off of the optical 
lattice is fast compared to the oscillation period, but slow 
enough to avoid diffraction of the atoms (i.e., adiabatic 
with respect to band excitations). 

We observe damped dipole oscillations for axial lattice 
depths from V — QEr to 2Er, as seen in Fig. ^ I n 
the absence of an axial lattice, we observe oscillations 
(period T — 15.4 ms) consistent with no damping (Fig. 
la), indicating that tube-to-tube dephasing and trap an- 
harmonicities are not significant on the timescale of our 
experiments. However, the oscillations are noticeably 
damped in a lattice only 0.25-Er deep. Such a shallow 
lattice modulates the atomic density by only 6%, and 
modifies the single-particle energy-quasimomentum dis- 
persion relation E(q) from that of a free-particle around 
only the last few percent of the BZ. (We note here, and 
discuss further below, that the amplitude of motion is 



kept well within the quadratic part of E(q) for shallow 
lattices). 

Beyond a lattice depth of « 3-Er the motion is over- 
damped, and there are no oscillations. In this case, the 
atoms' velocity can be quite small, so we use a tech- 
nique that maps the atoms' position in the trap to the 
cloud position after TOF. The experiment proceeds as 
before, except that after the trap is displaced by Zo, the 
atoms are allowed to relax toward their equilibrium po- 
sition at z — for a fixed time t w = 90 ms. We then 
rapidly (with time constant ~ 250 fis) turn off only the 
axial lattice. The remaining transverse lattice and mag- 
netic potentials are left on for 3.75 ms (approximately 
a quarter-period of undamped axial harmonic motion), 
then turned off simultaneously (as in the underdamped 
experiment). This converts the axial displacement z(t w ) 
into a velocity, which we measure by TOF. 

Figure 0i shows z(t w = 90 ms) as a function of axial 
lattice depth. For the shallowest lattices this wait time is 
sufficient for the atoms to damp to the equilibrium posi- 
tion z — 0. For the deepest lattices the motion is so over- 
damped that there is negligible motion during this time, 
and the position remains z w zq. We note that in the 
absence of damping, atoms would tunnel through the lat- 
tice to the equilibrium position in a time (T/4)y'm*/m, 
where m* is the effective mass. This time is only 8 ms 
for non-interacting particles in a IOEr lattice. 

To quantify the damping we model the motion as 
damped simple harmonic, m*z = —bz — kz, and extract 
a damping constant b = b(V) for different axial lattice 
depths V. For underdamped motion, we simultaneously 
fit the oscillation data for 8 depths to the expression 



z{t) = A- 



-bt/2m" 



sin(o^), 



where ui = yjk/m* — (b/2m*) 2 , and A and k are fit pa- 
rameters common across all V. For overdamped motion, 
we determine b = b(V) from the overdamped solution for 
z(t w ), which in the limit of strong damping simplifies to 



-kt„/b 



-bt m /r. 



1 — km* /b 2 1 — b 2 /km* 



where z and k are inputs derived from measurements of 
undamped oscillations. In our analysis we use a single- 
particle calculation of the effective mass m* [2^] . 

Figure |3 shows a plot of b(V)/bo versus lattice depth 
V, where 60 = 2mwo corresponds to critically damped 
harmonic motion for u>q = ^Jk/m. We show data for 
both the under- and overdamped regimes, and note that 
the damping constant increases by at least a factor 1000 
for a 30-fold increase in lattice depth. 

The axial width of the cloud after TOF can provide 
information about the distribution of atomic quasimo- 
menta in the lattice. The lattice turn-off time constant 
of 250 [is is long enough to avoid diffraction, but short 
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FIG. 2: Overdamped motion of a ID Bose gas in an optical 
lattice, (a) Plot of the atoms' position 90 ms after shifting the 
trap, as a function of lattice depth. Inset depicts relaxation of 
the atoms toward equilibrium. Immediately after the trap is 
displaced to zo, the atoms (open circle) begin to move toward 
equilibrium, reaching a displacement z from equilibrium after 
90 ms (solid circle). Also shown is the initial position of the 
atoms (dash-dot line), (b) The 1/e half- widths of Gaussian 
fits to axial TOF distributions (closed symbols), and the half- 
widths of square transverse TOF distributions (open symbols) 
resulting from a uniformly filled BZ (see text, and Fig. 2J. 
Also shown (dashed line) is the BZ calculated from lattice 
parameters. 



enough to be non-adiabatic with respect to inter- well tun- 
neling and interactions. (Related experiments in 2D [2lj 
and 3D |2y| lattices support this conclusion.) In the ab- 
sence of interactions, and neglecting the initial size of the 
cloud, the turn-off maps the single-particle quasimomen- 
tum distribution of atoms in the lattice to free particle 
(i.e., plane- wave) momentum states that can be directly 
observed in TOF. In the presence of interactions, the 
mapping is complicated by mean-field repulsion during 
TOF. A variational calculation [53 indicates that mean- 
field repulsion is in fact the dominant contributor to the 
axial TOF width in our system. Therefore, the extracted 
TOF width greatly overestimates the width of a narrow 
initial quasimomentum distribution. 

An example TOF image is shown in Fig. together 
with cross-sectional profiles of the optical depth along 
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FIG. 3: Plot of the reduced damping constant b/bo for 
various depths V of the axial lattice, as determined using 
the underdamped (squares) and overdamped (circles) exper- 
imental techniques. Near critical damping (V ~ 3Er), the 
analysis cannot distinguish between underdamped and over- 
damped motion, so only upper and lower bounds are shown 
at V = 2.25 and 2.50. 
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FIG. 4: Cross-sections (a) of a TOF absorption image (b) of 
the expanded atom cloud. The plot of the transverse cross- 
section (along x) is offset vertically by 0.5 units for clarity. 
The solid lines are Gaussian (triangular) fits to the axial 
(transverse) profiles, respectively. The dashed ovals in (b) in- 
dicate the peak-to-peak range of motion of undamped dipole 
oscillations. 



the axial and transverse directions. The first BZ of the 
transverse lattice is uniformly filled, producing a uniform, 
square spatial distribution (in the xy plane) after TOF. 
Our imaging system views this square distribution along 
the diagonal in the xy plane, resulting in a triangular 
profile, from which we extract the width of the square 
(open symbols, Fig. [2t>) ■ In the axial direction the distri- 
bution is narrower and reasonably well fit to a Gaussian, 
from which we extract the axial (along z) TOF 1/e half- 
width (filled symbols, Fig. \Bp). Even for the strongly 
overdamped data, the axial TOF width (which, we recall, 
overstates the width of narrow quasimomentum distribu- 
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tions), is much narrower than the BZ. This implies that 
the inhibition of transport is not due to effects related 
to Bloch oscillations of a filled BZ, as observed in Refs. 

Furthermore, we do not see a significant difference 
in TOF width between atoms that undergo damped har- 
monic motion and those that are unexcited but held for 
an equal time 28]. This is in stark contrast to earlier 
experiments on 3D BECs 0, @, Q] , where strong damp- 
ing was accompanied by a pronounced broadening and 
fragmentation of the quasimomentum distribution. 

Large amplitude dipole oscillations in a lattice can 
damp due to dynamical instabilities caused by particle 
interactions. For a 3D BEC moving in an optical lattice, 
such an instability point occurs at q > q^/2 (q^ = 2ir/X 
is the BZ boundary), where the dispersion relation has an 
inflection point, as predicted in [SUlfJ, and observed in 
IE HL E3 • This effect is manifested as a large increase in 
the width of the quasimomentum distribution. Here, in 
contrast, we keep the maximum (single-particle) quasi- 
momentum g max of the oscillation small by limiting the 
initial energy of displacement E(zq) — itiu!qZq/2. For 
V < 2Er, our choice of zq corresponds to g max ~ <Ztt/5. 
For deeper lattices our fixed zq corresponds to a larger 
<Zmax, but is always less than q n /2 for V < 9Er. (In a 
separate experiment, we excited oscillations in our sys- 
tem with twice the usual amplitude, and saw stronger 
damping that was accompanied by a broadening of the 
axial TOF width by nearly a factor of two.) 

For small amplitude dipole oscillations (g max <S q-n/2) 
of a 3D BEC in an optical lattice, the effect of the lat- 
tice is merely to increase the effective mass, leading to 
undamped motion at a lower frequency HQ. In the 
reduced dimensionality system of our ID Bose gas, we 
have seen that the optical lattice has a qualitatively dif- 
ferent effect. To highlight the difference between these 
two situations, we excited dipole oscillations in a 3D BEC 
(i.e., no transverse confining lattice) in a 4Er axial lat- 
tice, and saw no damping. This is in contrast to the re- 
sults of the same experiment in a ID Bose gas (i.e., with 
a transverse confining lattice), shown in Fig. |3J where 
b(V = 4£tj)/&o « 50 corresponds to extremely over- 
damped motion. 

After we performed these experiments, theoreti- 
cal treatments appeared which suggested that zero- 
temperature quantum fluctuations can lead to substan- 
tial damping of transport in a ID atomic Bose gas 
[2H Isflj . Our observations, including the significant 
damping in lattices too shallow to support a Mott- 
Insulator (MI) phase l31l , can be explained by the mech- 
anisms of Refs. [2|| |3(|, but appear to be inconsistent 
with a mechanism involving incompressibility |l3| . 

It is possible that there is a temperature dependence 
to the damping; unfortunately, we can derive little in- 
formation on temperature from the TOF widths. Future 
experiments to investigate the temperature dependence 
could shed light on the relative importance of quantum 



and thermal fluctuations to dissipation, a question of in- 
terest in, for example, the development of ultrathin su- 
perconducting wires |l7| . We also look forward to test- 
ing other explicit predictions of these theories, such as 
the dependence of damping on displacement and dimen- 
sionality. We note that the periodic potential of an opti- 
cal lattice is free from defects, and that ID atomic Bose 
gases are well isolated from the environment, yielding a 
relatively clean system in which to compare experiment 
with theory. The ability to continuously and dynami- 
cally vary the confining potentials makes optical lattice 
experiments attractive for future studies of superfluidity 
in low dimensional quantum systems. 

The authors gratefully acknowledge helpful discussions 
with Ana-Maria Rey, Julio Gea-Banacloche, and Guido 
Pupillo. This work was partially supported by ARDA. 



* Current address: Department of Physics, Pennsylvania 
State University, University Park, PA 16802 

[1] G. Modugno, F. Ferlaino, R. Heidemann, G. Roati, and 
M. Inguscio, Phys. Rev. A 68, 011601(R) (2003). 

[2] L. Pezze, L. Pitaevskii, A. Smerzi, S. Stringari, G. Mod- 
ugno, E. DeMirandes, F. Ferlaino, H. Ott, G. Roati, and 
M. Inguscio, Phys. Rev. Lett. 93, 120401 (2004). 

[3] H. Ott, E. de Mirandes, F. Ferlaino, G. Roati, G. Mod- 
ugno, and M. Inguscio, Phys. Rev. Lett. 92, 160601 
(2004). 

[4] O. Morsch, J. H. Miiller, M. Cristiani, D. Ciampini, and 
E. Arimondo, Phys. Rev. Lett. 87, 140402 (2001). 

[5] F. S. Cataliotti, S. Burger, C. Fort, P. Madaloni, F. Mi- 
nardi, A. Trombettoni, A. Smerzi, and M. Inguscio, Sci- 
ence 293, 843 (2001). 

[6] S. Burger, F. S. Cataliotti, C. Fort, F. Minardi, M. In- 
guscio, M. L. Chiofalo, and M. P. Tosi, Phys. Rev. Lett. 
86, 4447 (2001), b. Wu and Q. Niu, Phys. Rev. Lett. 89, 
088901 (2002), and S. Burger, et al., ibid., 088902 (2002). 

[7] F. Ferlaino, P. Maddaloni, S. Burger, F. S. Cataliotti, 
C. Fort, M. Modugno, and M. Inguscio, Phys. Rev. A 
66, 011604(R) (2002). 

[8] F. S. Cataliotti, L. Fallani, F. Ferlaino, C. Fort, P. Mad- 
daloni, and M. Inguscio, New J. Phys. 5, 71 (2003). 

[9] B. Wu and Q. Niu, Phys. Rev. A 64, 061603(R) (2001). 
[10] A. Smerzi, A. Trombettoni, P. G. Kevrekidis, and A. R. 

Bishop, Phys. Rev. Lett. 89, 170402 (2002). 
[11] L. Fallani, L. DeSarlo, J. E. Lye, M. Modugno, R. Saers, 
C. Fort, and M. Inguscio, Phys. Rev. Lett. 93, 140406 
(2004). 

[12] B. L. Tolra, K. M. O'Hara, J. H. Huckans, W. D. Phillips, 
S. L. Rolston, and J. V. Porto, Phys. Rev. Lett. 92, 
190401 (2004). 

[13] T. Stoferle, H. Moritz, C. Schori, M. Kohl, and 

T. Esslinger, Phys. Rev. Lett. 92, 130403 (2004). 
[14] H. Moritz, T. Stoferle, M. Kohl, and T. Esslinger, Phys. 

Rev. Lett. 91, 250402 (2003). 
[15] B. Parades, A. Widera, V. Murg, O. Mandel, S. Foiling, 

I. Cirac, G. V. Shylapnikov, T. W. Hansen, and I. Bloch, 

Nature 429, 277 (2004). 
[16] T. Kinoshita, T. Wenger, and D. S. Weiss, Science 305, 



5 



1125 (2004). 

[17] A. Bezryadin, C. N. Lau, and M. Tinkham, Nature 404, 
971 (2000). 

[18] F. Dalfovo, S. Giorgini, L. Pitaevskii, and S. Stringari, 
Rev. Mod. Phys. 71, 463 (1999). 

[19] M. Kramer, L. Pitaevskii, and S. Stringari, Phys. Rev. 
Lett. 88, 180404 (2002). 

[20] M. BenDahan, E. Peik, J. Reichel, Y. Castin, and C. Sa- 
lomon, Phys. Rev. Lett. 76, 4508 (1996). 

[21] M. Greiner, I. Bloch, O. Mandel, T. W. Hansch, and 
T. Esslinger, Phys. Rev. Lett. 87, 160405 (2003). 

[22] Error bars here and everywhere represent the uncorre- 
lated combination of 1 standard deviation statistical and 
systematic uncertainties (if any). 

[23] S. L. Winoto, M. T. DePue, N. E. Bramall, and D. S. 
Weiss, Phys. Rev. A 59, R19 (1999). 

[24] For all lattice depths, the maximum site-to-site potential 
difference of 40 Hz is much less than both the on-site 
interaction energy and the tunneling energy, so that the 
transport studied in this work is fundamentally different 
from the excitations due to potential gradients reported 
in Ref. .26]. 



[25] For overdamped motion, the extracted damping con- 
stants b depend very weakly on m* /m. For underdamped 
motion, the shallowness of the axial lattice means correc- 
tions to m* due to interactions should be small. There- 
fore we expect all extracted damping constants b to be 
insensitive to corrections to m* due to interactions. 

[26] M. Greiner, O. Mandel, T. Esslinger, T. W. Hansch, and 
I. Bloch, Nature 415, 39 (2002). 

[27] A. M. Rey, private communication. 

[28] During the ~ 100 ms in which we observe the cloud's mo- 
tion, we observe a small growth in width of the TOF im- 
age, at a rate that is independent of excitation or damp- 
ing, consistent with heating due to spontaneous scatter- 
ing of lattice photons. 

[29] A. Polkovnikov and D.-W. Wang, Phys. Rev. Lett. 93, 
070401 (2004). 

[30] J. Gea-Banacloche, A. M. Rey, G. Pupillo, C. J. 

Williams, and C. W. Clark (2004), cond-mat/0410677. 
[31] C. Kollath, U. Schollwock, J. von Delft, and W. Zwerger, 

Phys. Rev. A 69, 031601 (R) (2004). 



